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Abstract 

We show that GLAST should be able to support or rule out the unresolved-blazar hypothesis for the 
origin of the diffuse, extragalactic gamma ray background. 

1 Introduction: 

Q"^ ' Perhaps the most promising model proposed for the origin of the ~GeV extragalactic 7-ray back- 

0^ ■ ground (EGRB), first detected by SAS-2 [Fichtel, Simpson, & Thompson 1978] and later confirmed 
' by EGRET, is that it is the collective emission of an isotropic distribution of faint, unresolved blazars 
[Stecker & Salamon 1996 and references therein]. EGRET has determined that the EGRB spectrum 
Q ■ is consistent with a single power-law, 
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\ between 0.1 and ~ 50 GeV (statistics limited) [Sreekumar et al. 1998]. Because GLAST has a 
' roughly energy-independent area of 10^ cm^ above 0.1 GeV (compared to EGRET's 10'^ cm^ at 1 
GeV and 10^ cm^ at 100 GeV), with an estimated point source sensitivity (PSS) nearly two orders 
\ of magnitude lower than EGRET's, GLAST will be able to (a) detect something on the order of 10^ 
Q.^ ' times more blazars than EGRET, and (b) measure the EGRB spectrum to > 1 TeV (assuming the 
0^ ■ EGRET power law spectrum). These two capabilities will enable GLAST to either strongly support 
or reject the unresolved-blazar hypothesis for the origin of the EGRB. 

^ \ 2 The Unresolved Blazar Model: 

\ To determine the collective output of all 7-ray blazars, one can use the observed EGRET distribu- 

^ ■ tion of 7-ray luminosities and extrapolate to obtain a "direct" 7-ray luminosity funtion (LF) f-y{l-y, z) 
[Chiang, et al. 1995 ] (/-y units being sources/co-moving volume-differential luminosity), where is 
the differential luminosity (s~^) at a (source) fiducial photon energy Ef, and we assume power law 
^ ' spectra for all sources, 1{E) = 1^{E/Ef)~°'. Alternatively, one can make use of much larger catalogs 
■ at other wavelengths, and assume some relationship between the source luminosities at the catalog 
wavelength and the GeV region [Padovani et al. 1993; Stecker, Salamon, &; Malkan 1993]. Both 
methods are fraught with potential errors. In the former method, only the "tip of the iceberg" of 
the 7-ray LF has been observed by EGRET, and extrapolation to lower luminosities that fall below 
EGRET's PSS (for any reasonable source distance) involves some level of assumption. In the latter 
method, the assumption of a linear relation between the luminosities of a source at disparate wave- 
lengths is by no means well established [Padovani, et al. 1993; Muecke, et al 1996, Mattox, et al. 
1997]. 

We used the latter method in the past [Stecker and Salamon 1996] to estimate the contribution of 
unresolved blazars to the EGRB, and found that up to 100% of EGRET's EGRB can be accounted 
for. This model assumes a linear relationship between the differential 7-ray luminosity at Ef = 0.1 
GeV and the differential radio luminosity Ir at 2.7 GHz for all sources, viz., = nlr, where k is a 
constant. One can then use the measured radio LF fr{lr,z) for blazars (primarily for flat spectrum 
radio quasars) [Dunlop &; Peacock 1990] to calculate the collective 7-ray output of all blazars. 



The simplified elements of the calculation follow (the details can be found in Stecker and Salamon 
1996). First, the 7-ray and radio LFs are related by f.y{l^,z) = fr{K~^l^, z). The number of 
sources TV detected by a detector is a function of the detector's point source sensitivity (PSS) at the 
fiducial energy Ef, [F{Ef)]tnin, where the integral 7-ray number flux F is related to l-y by 

F{E) = l^{E/Efy/ATTa{l + a)''+^Rlr^, (2) 

where R^r is the luminosity distance to the source. The number of sources at redshift z seen at Earth 
with an integral flux F(Ef) is given by 



dF{Ef) 



AF{Ef) = J ATTRlr'^dr f^{ly,z{r))Al^, (3) 



where l-y in the integrand depends on z{r) and F{Ef) from Eq.§. Figure 1 shows our 1996 calcu- 
lations for the number of sources versus flux, compared to the EGRET detections. The cutoff at 
10~'^cm~^s~^ for Ef = 0.1 GeV, their quoted PSS, is evident. 
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Figure 1: Source number count per one- fifth decade of integral flux at Earth. The straight dotted 
line is the Euclidean relation Af{> F) oc F~^/'^ for homogeneous distribution of sources. The open 
circles represent the EGRET blazar detections, and the solid line is the model prediction. 



To calculate the EGRB, we integrate over all sources not detectable by the telescope to obtain the 
differential number flux of EGRB photons at an observed energy Eq: 

dN^ , „ , /",„•:! 9 , /■ , .X Z"'™"^ dF 
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iTTRy dr J dapia) ^(^o(l + z))!^^, z)e-^(^0'^) dL,. (4) 



This expression includes an integration over the probabihty distribution of spectral indices a (based 
on the 2nd EGRET Catalog [Thompson et al. 1995]). There is also an important attenuation factor 
in this expression, due to the loss of 7-rays as they propagate through the intcrgalactic medium and 
interact with cosmic UV, optical, and IR background photons to produce pairs. Although there 
are significant uncertainties in estimates of the z-dependence of the soft photon background [Salamon 
&; Stecker 1998], a qualitative feature cannot be avoided: If a substantial fraction of the EGRB is 
from high- 2: flat spectrum radio quasars, a steepening in the spectrTim should be seen at energies 
above 20 GeV. Figure 2 shows the calculated EGRB spectrum (based on EGRET's PSS) compared 
to EGRET data. The "bump" in the spectrum above 10 GeV is due to there being a finite width in 
the assumed power-law spectral index distribution of blazars, which produces a summed unabsorbed 
spectrum with a positive second derivative [Stecker & Salamon 1996]. Although the match between 
the model's spectrum and the recent EGRET EGRB data [Srcckumar et al. 1998] does not appear 
to be very good, one must examine it the light of the assumptions and uncertainties. The location 
and amplitude of the "bump" depend upon the amount of high redshift absorption, which, in turn, 
depend upon uncertain UV and optical soft photon background densities [Salamon & Stecker 1998]. 
The "bump" also depends upon the poorly known blazar spectral index distribution, which we have 
assumed to independent of quiescent-state luminosity. Finally, we not that the assumed power-law 
form for high redshift quasars above 10 GeV energy is highly questionable in light of the predictions 
of the Compton models for 7-ray emission from these sources; in fact their spectra are predicted 
to steepen at these energies. We also note that EGRET's EGRB error bars are systematic, not 
statistical and that Sreekumar et al. [1998] have stated that the model cannot be ruled out based on 
the EGRET data. 

3 GLAST and the EGRB: 

Figure 1 shows that O(IO^) blazars should be detectable by GLAST, assuming it achieves a PSS 
of ~ 2 X 10"^ cm~^s~^, and Figure 2 shows that this should reduce the EGRB by a factor of ~ 2 for 
E>1 GeV. These predictions are arrived at by replacing EGRET's PSS with GLAST's in our 1996 
model. 

Figure 2 shows our predicted reduction in the EGRB from unresoved blazars to be observed by 
GLAST compared with that observed by EGRET. The dashed line in Figure 2 represents the energy 
interval in which it is difficult to determine GLAST's EGRB spectrum, due to the eff'ect of worsening 
angular resolution with decreasing energy. For a lower PSS level, the number of detectable sources 
is larger. However, the angular resolution limits the number of distinct sources one can identify, 
crudely given hy Af{E) 4Tr/Trag{E), where (Tg{E) is the energy-dependent angular resolution. When 
^(> Fmin) > J^{E), individual sources are no longer separable from the EGRB. With GLAST's 
proposed angular resolution function (GLAST Science Document), the equality 7V(> Fmin) = M{E) 
occurs at E somwhat less than 1 GeV; thus ioi E < 1 GeV there is not as much a reduction in the 
level of the EGRB (compared to EGRET) as for energies E >l GeV. 

We conclude that GLAST can test the unresolved blazar background model in three ways: (1) 
GLAST should sec roughly 2 orders of magnitude more blazars than EGRET because of its ability 
to detect the fainter blazars which contribute to the EGRB in our model. Moreover, GLAST can 
test our assumption of an average linear relation between 7-ray and radio flux and our assumed 
redshift distribution of blazars by testing the details of our source count versus flux prediction. (2) 
With GLAST's improved PSS leading to more blazars being resolved out, fewer unresolved blazars 
will be left to contribute to the EGRB, thus reducing the level of the measured EGRB compared to 
EGRET's. (3) GLAST's much greater aperture at 100 GeV will allow a determination of whether or 
not a steepening exists in the EGRB, since the number of EGRB 7-rays recorded by GLAST above 
100 GeV will be of order 10^ to 10^, assuming a continuation of the EGRET power law spectrum. 
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Figure 2: The upper and lower solid curves correspond to model predictions for the EGRB spectrum 
given EGRET'S and GLAST's PSS (respectively). The data points are from the latest EGRET 
EGRB analysis [Sreekumar ct al. 1998]. The dashed line segment in the lower curve corresponds to 
the energy interval where increasing source confusion, owing to worsening angular resolution at lower 
energies, prevents the discernment of point sources, which then contribute to the EGRB. 



However, this last test must be qualified because of the unknown amount of steepening in the high 
reshift quasar spectra above 10 GeV which can mimic an absorption effect. 
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